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ABSTRACT 

We examine the recent results of the MACHO collaboration towards the 
Large Magellanic Cloud (Alcock et al. 1996) in terms of a halo brown dwarf or 
white dwarf population. The possibility for most of the microlensing events to 
be due to brown dwarfs is totally excluded by large-scale kinematic properties. 
The white dwarf scenario is examined in details in the context of the most recent 
white dwarf cooling theory (Segretain et al. 1994) which includes explicitely the 
extra source of energy due to carbon-oxygen differentiation at crystallization, 
and the subsequent Debye cooling. We show that the observational constraints 
arising from the luminosity function of high-velocity white dwarfs in the solar 
neighborhood and from the recent HST deep field counts are consistent with a 
white dwarf contribution to the halo missing mass as large as 50%, provided i) 
an IMF strongly peaked around ~ 1.7 M and ii) a halo age older than ~ 18 
Gyr. 

Subject headings : stars: low-mass, brown dwarfs — stars: white dwarfs - 
stars : luminosity function, mass function — The Galaxy : halo — dark matter 
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1. Introduction 

There is compelling evidence for believing that a large amount, if not most of baryonic 
matter in the Universe is under the form of dark, unobserved objects. On the other hand, 
there is also evidence that spiral galaxies are surrounded by a large amount of non-luminous 
mass which is responsible for their observed small-scale and large-scale kinematic properties 
(velocity dispersion and rotation curve). These two facts yield the natural conclusion that 
baryonic dark matter is a plausible candidate for halo dark matter. A breakthrough in this 
longstanding, unsolved problem has been accomplished recently with the developement of 
microlensing experiments, by inferring the presence of dark objects in the halo through their 
gravitational effect on luminous matter. A detailed analysis of the first year of the MACHO 
(Alcock et al. 1993) and EROS (Aubourg et al. 1993) observations towards the LMC, 
complemented by the determination of the mass function of low-mass halo stars, yielded an 
average mass for the dark objects < m >^ 0.03 M , well within the brown dwarf domain 
(Mera, Chabrier and Schaeffer 1996a). The inferred maximum brown dwarf contribution to 
the halo mass budget was found to be ~ 10 to 20% (Alcock et al. 1995; Gates, Gyuk and 
Turner, 1995; Mera et al. 1996a). These results have now to be reconsidered in the light of 
the most recent analysis of the MACHO collaboration, which includes one more year LMC 
data. This yields now a total of seven candidates for MACHO with longer durations, from 
30 to 110 days (Alcock et al. 1996). 

In this Letter, we re-examine the M-dwarf, brown dwarf and white dwarf contributions 
to the halo mass budget in light of the most recent observational constraints (star counts 
and MACHO results) and white dwarf cooling theory. 
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2. M-dwarfs, brown dwarfs 

HST star counts at large magnitudes (Bahcall et al. 1994, Hu et al. 1994) show that 
the M-dwarf contribution to the Galactic missing mass M^ yn « 10 12 M & represents at 
most a few percents. A more precise determination can be obtained from the observed 
luminosity function (LF) of high-velocity M-dwarfs in the solar neighborhood (Dahn et al. 
1995). A detailed analysis of this LF shows that the mass function (MF) of the spheroid 
(characterized by a ~ 1/r 3 density-profile) is reasonably well described by a power-law 
MF 4>(m) oc m~ a with a ~ 2 — 2.5, from 0.6 M & down to the hydrogen burning limit 
and a normalization dN/dm(0.1 M ) = 1CT 2 ' 7 M^pc -3 (Mera et al. 1996a). This yields 
a maximum local density for the spheroid+halo M-dwarf population ~4x 1O _5 M pc~ 3 , 
and thus an optical depth TMdwarf < 4 x 10~ 9 . A comparison with the value inferred 
from the new MACHO results r « 1.7 X 1(T 7 shows convincingly that M-dwarfs can be 
responsible for less than 0.1 of the microlensing events towards the LMC On the other 
hand, reasonnable estimates for the LMC characteristics yield an optical depth r«5x 10~ 8 , 
a factor 3 to 4 smaller than the recent MACHO results (Sahu 1994). 

As mentioned above, the analysis of the previous MACHO+EROS results showed that 
the observed events were likely to be due to halo brown dwarfs, with an average mass 

< m >m 0.03 M (Mera et al. 1996a). These calculations must be re-examined in the 
context of the recent MACHO results. The new 7 events yield < t e >oc< ^pm >< - 40 
days (Alcock et al. 1996)[]. Since the minimal tangential velocity < v > of the lens is bound 
by the rotation velocity of the line of sight (~ 220 km s _1 )(Mera et al. 1996b), this excludes 
totally the possibility for these events to be due to brown dwarfs and yields an average mass 

< m >m 0.5 M . Since M-dwarfs are excluded, the inferred optical depth r w 1.7 x 10~ 7 



*Note that the event time-scale is defined here as the Einstein radius crossing time, 
whereas the MACHO group adopts the Einstein diameter crossing time in their definitions. 
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means that about 40% of the hidden mass consists of halo white dwarfs. [] 

3. White dwarfs 

The white dwarf contribution to the halo missing mass has been examined over the 
past by different authors. These studies focussed mainly on the constraint arising from the 
Galaxy chemical evolution (Olive 1986; Ryu et al. 1990). These authors show that WDs 
are unlikely candidates for providing the entire halo missing mass. Recent calculation of the 
expected radiation signature of the progenitors in galactic halos at large redshift (Chariot 
and Silk 1995) show that a white dwarf mass fraction larger than ~ 10% of the missing 
mass would be in conflict with observations. These results will be commented below. 

An independent, more stringent constraint, comes from the observed white dwarf 
luminosity function (WDLF) in the solar neighborhood, as considered by Tamanaha, 
Silk, Wood and Winget (1990). However, although pointing the way, these calculations 
were based on simplified WD interior and a WD cooling theory aimed at describing the 
disk WDLF, thus appropriate for objects younger, and thus warmer, than the expected 
halo population. In particular these calculations do not include a complete treatment of 
crystallization (see below), which occurs around logL/L Q rs —3.5 in WD interiors. This 
affects substantially the cooling of halo WDs and will modify significantly the expected 
halo WDLF. Q 

^Neutron stars and stellar black holes must be rejected as a significant halo population 
mainly on the basis of the severe constraints arising from the observed metallicity and helium 
abundances (Ryu, Olive and Silk 1990). In any case neutron stars and black holes are likely 
to contribute considerably less than WDs to the dark mass fraction. 

* After the present calculations were completed, we were aware of similar studies by Adams 
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In this Letter, we use the most updated WD cooling theory for carbon/oxygen WDs, 
with the appropriate equation of state both in the classical and in the quantum (crystal) 
regime (Segretain et al., 1994; Chabrier, 1993) and a helium-rich atmosphere (Wood 1992), 
characteristic of most cool (T cff < 6000 K) so-called "DB" WDs. As first suggested by 
Stevenson (1980), the gravitational energy release due to carbon-oxygen differentiation at 
crystallization affects drastically the subsequent cooling time of the star, thus changing 
the luminosity for a given age (Segretain & Chabrier 1993). A consistent treatment of the 
crystallization phase diagram along WD evolution has been derived recently by Segretain 
et al. (1994). As shown by these authors, the crystallization processes modify appreciably 
the WD cooling time and then the WDLF for logL/L ^ —4, characteristic of old disk and 
halo WDs. The LF derived with this theory yields an estimate for the age of the Galactic 
disk Tdisk ~ 10.5 — 12 Gyr, depending on the bolometric correction used for the observed LF 
(Hernanz et al. 1994), about 20% larger than estimates based on cooling theories which do 
not include the complete crystallization process (Wood 1992; see Segretain et al. 1994 §4.1 
for details). 

The calculations proceed as in Hernanz et al. (1994). The WDLF reads : 



Here t coo \ = dt coo i/ 'dM^ is the characteristic cooling time, where t coo i is the WD 
cooling time. t ms and th denote respectively the age spent on the main sequence for the WD 
progenitor and the age of the halo. The function 0(m) is the initial mass function and m in f 
and m sup denote respectively the minimum and the maximum mass of the WD progenitors 

and Laughlin (1996). These calculations, however, are based on crude (pure carbon) WD 
interior and the afore-mentioned cooling theory. 




r cool (L,m) x ip[t h - t cool (L,m) - t ms (m)] x <p(m)dm 



(1) 
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which contribute at luminosity L. Since the age of the halo is much larger than any time 
associated with star formation, the initial stellar formation rate ip(t) is well approximated 
by a burst at t — 0, i.e. a 8(t = 0) function. In that case eqn(l) reduces to : 

where v(th — t coo i) represents the number of WDs formed at t = th — t coo i, i.e. the 
number of stars with a main sequence lifetime t ms = th — t coo i. We verified that finite-time 
SFR, e.g. a constant SFR along At ^ (a reasonnable representation of a continous series 
of burst SFRs) or an exponential SFR yield very similar results. The progenitor- WD mass 
relation is m WD 0.45 + 0.1m (Iben and Tutukov 1984). The WDLF is normalized to : 



J ndM bo i = -2.5 J n dlog(L/L ) = X WD p dyn / < m WD > pc 3 (3) 

where X WD is the (sought) mass fraction under the form of WDs in the halo of the 
Galaxy. As shown in eqn(2) the most essential parameter in this calculation is the white 
dwarf cooling time t coo i. We use the afore-mentioned WD cooling sequences calculated in 
Segretain et al. (1994) and Garcia-Berro et al. (1996). 

The second important parameter to be determined in Eq.(l) is the IMF 0(m). As 
mentioned in §2, a severe constraint arises from the recently determined mass-function 
(slope and normalization) of halo M-dwarfs. The predicted star counts obtained with 
this MF for a spheroid(l/r 3 )+halo(l/r 2 ) density profile are in perfect agreement with 
the observations of the HST at large magnitude (I > 25) (Mera et al., 1996a). On the 
other hand, the observed halo metallicity implies that stars above m ^ 8 M , believed 
to be type II Supernovae progenitors, represent at most ~ 1% of the halo initial stellar 
population (Ryu et al. 1990). These observational constraints show that, for WDs to 
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contribute significantly to the mass of the halo, the IMF must exhibit a strongly bimodal 
behaviour and peak around some characteristic mass in the range [m^/, ~ 8M ]. The 
minimum mass corresponds to a main-sequence lifetime of the progenitor equal to the age 
of the halo, i.e. m in f ps 0.9 M for t — 10 to 25 Gyr. Several functional forms can be 
advocated for such an (unknown) IMF. We elected a simple cut-off power-law function 
4>(m) = dN/dm = A e~( m l m ) (ix mT^' 1 (see e.g. Larson 1986). This form mimics adequatly 
a strongly peaked IMF and is very similar to functional forms based on stellar formation 
theory (Adams and Laughlin 1996). The IMF is normalized to : 



which determines A (for a given Xwd and mwDijn) relation). The parameter-space for 
(m, /3i, fa) is constrained by the required negligible number of stars outside the mass-range 
[~ O.9M ,~ 8M ] but different values yield quantitatively different mass-distributions. A 
large number of masses > 2 M would raise severe problems for the fraction of ejected gas 
and the subsequent helium and metal galactic enrichment (Hegyi and Olive 1986; Ryu et al. 
1990). In order to examine the dependence of the IMF on the results, we thus considered 
two functions, namely (m = 2.0, /?i = 2.2, /3 2 =5.15), peaked around ~ 1.3 M (hereafter 
IMF1), and (m = 2.7, fa = 2.2, #2=5.75), peaked around ~1.7M (hereafter IMF2). The 
complete M-dwarf+WD IMF's fulfilling all the afore-mentioned constraints are shown on 



Observational constraints. The LF of field WDs has been obtained by Liebert et al. 
(1988) up to My 19 (i.e. L/L Q <; 10~ 5 — 10~ 6 , depending on the bolometric correction 
BCy). The LF declines abruptdely for My ps 16, which corresponds to logL/L ps —4.2 
to —4.6. As stated by these authors, no WD was found at fainter magnitudes, with this 
or with other proper-motion samples, whereas stars up to My = 19, i.e. three magnitudes 




Figure 1. 
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fainter, have been observed with similar programs (Liebert et al. 1988; Monet et al. 1992). 
Interestingly enough, five WDs in the Liebert et al. sample have tangential velocities 
Vtan > 250 kms^ 1 and My > 13 and thus are assignable to the halo sample (shown by filled 
circles on Figure 2). 

More recent constraints arise from the HST counts up to I = 26.3 (Flynn et al. 1996). 
For a WD mass fraction Xwd ~ 10%, it is easy to show that the observed number of WDs 
in the HST field implies that halo WDs must have My « M boi ^ 14. This is consistent with 
the observed high-velocity WDs. The colors of the disk and halo WDs were taken to be 
^ V — I ^ 2, as suggested by the observations (Liebert et al. 1988; vonHippel, Gilmore 
& Jones 1995), and the bolometric correction ^ BCy ^ 1 (Liebert et al. 1988; Bergeron, 
Saumon & Wesemael 1995). 

The observed WDLF is represented on Figure 2, with different halo WDLFs. The 
dotted line is the disk WDLF from Segretain et al. (1994) for an age = 10.5 Gyr 0. The 
crosses correspond to a 90% exclusion confidence level in the limit of detection, i.e. the 
possibility to see at least two WDs above this line whereas none has been detected is rejected 
at the 2-<j level (see e.g. §IV of Liebert et al. 1988). The solid lines show the halo WDLFs 
for halo ages t h = 14, 16, 18 and 20 Gyr, normalized to X WD = 1, 2, 4 and 8% respectively, 
for calculations done with IMF1. For the distribution of progenitors corresponding to this 
IMF, differentiation at crystallization in the WD interiors leads to a bump in the halo 
WDLF in the range —5 ^ log(L/L Q ) ^ —4, therefore ruling out substantial WD mass 
fractions. This shows convincingly the importance of a complete treatment of crystallization 
in WD cooling. Calculations with no carbon/oxygen differentiation will underestimate the 

^based on a blackbody bolometric correction BCy for the observed LF (see Liebert et 
al., 1988)). A zero BCy will yield ~ 1.5 Gyr older ages (see Hernanz et al. 1994), for both 
the disk and the halo. 
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number of WDs by more than a factor ~ 5, for a given age and luminosity. Conversely 
they will yield halo ages ~ 2 Gyr younger for a given LF. In the same vein, an incorrect 
Debye treatment will change significantly the shape of the WDLF. The dashed lines 
correspond to the same calculations when using IMF2. The normalizations correspond now 
to X WD = 1.7,8,25 and 50%, for the same halo ages. Clearly, for the IMF2, a halo WD 
mass fraction <: 30%, in agreement with the MACHO results, can not be excluded, provided 
a halo age ^18 Gyr. 

We have compared the star counts predicted by these WDLF's with the recent HST 
observations at large magnitudes (Flynn et al. 1996), for a l/r 3 -spheroid and a 1/r 2 halo. 
All WDLFs predict at most (depending on BC V ) ~ 1.4 WD in the HST field at the limit 
magnitude / = 26 for a 100% WD halo, and thus are consistent with the HST counts. 
This shows that for these scarce and faint objects, large, nearby surveys put more severe 
constraints than deep pencil searches. 



4. Conclusion 

In this Letter, we have examined the possibility for the recent MACHO events 
(assuming these events are genuine microlensing events) to be due either to brown dwarfs or 
to white dwarfs. This determines directly the contribution of these objects to the missing 
mass in the halo of the Galaxy. Brown dwarfs are clearly excluded as a significant halo 
population. The luminosity function of halo white dwarfs has been calculated with the most 
accurate white dwarf cooling theory presently available. This WDLF is confronted to all 
available observational constraints on halo objects. We show that, under the two necessary 
conditions that i) the IMF in the halo differs totally from the one in the disk and exhibits a 
strongly peaked behaviour around m ~ 1.5 — 2 M , and ii) the halo is older than ~ 18 Gyr, 
the white dwarf mass fraction in the halo can represent ~ 25 to 50% of the dark matter 
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density, in agreement with the recent MACHO results. This would imply an initial stellar 
mass fraction > 50% and thus an essentially baryonic halo. These results are consistent 
with the ones obtained from galactic chemical evolution (Ryu et al. 1990), though they are 
in conflict with the conclusion raised by these authors that the disk must form no later than 
the halo. However, as stated by these authors, alternative scenarios in the disk formation 
can be advocated : the left-over gas fraction might have been ejected into the intergalactic 
medium, as suggested by recent observations of metal-rich hot gas in the Local Group (Suto 
et al. 1996)Q. The present results are also consistent with the ones obtained by Chariot and 
Silk (1995), based on the expected radiation signature in high-redshift galactic halos. These 
authors considered a Hubble time < 13 Gyr, and solar metallicity (i.e. slowly evolving) 
stars. The evolution of significantly older, i.e. highly redshifted, low-metallicity stellar 
populations will certainly be consistent with these observational constraints. 

Therefore, although providing a plausible explanation for the MACHO observations 
and the halo missing mass, the present scenario relies on the necessity to invoke a very 
peculiar, fine-tuned IMF. These calculations illustrate the difficulty to reconcile the recent 
MACHO results with other observational constraints. A detailed analysis of the OGLE, 
MACHO and EROS results, in the context of a consistent model for the Galaxy, will be 
presented in a forthcoming paper (Mera, Chabrier and Schaeffer 1996b). 



^Note also that, given the low-metallicity and thus the probably less efficient mass- loss, 
halo WDs may have larger masses than disk WDs, thus resulting in a smaller fraction of gas 
eject a 
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FIGURE CAPTIONS 

Figure 1 : Halo initial mass function. The dotted line illustrates the M-dwarf MF 
dM/dm oc m~ 2 - 2 (Mera et al. 1996a). The solid line is the IMF 0(m) = Ae" ( " i/mm J' )?1 m" ft 
with m = 2.0, ft = 2.2, #2=5.15 (IMF1). The dot-dashed line is the same IMF with 
m = 2.7, ft = 2.2, ft = 5.75 (IMF2). 

Figure 2 : White dwarf luminosity function (pc -3 M^). Empty circles : Liebert et al. 
(1998). Filled circles : high- velocity WDs (Liebert et al. 1988). Crosses : limit of detection 
at the 2-cr level (see text). Dotted line : disk WDLF for t d = 10.5 Gyr. Solid lines : halo 
WDLF for t h = 14, 16, 18,20 Gyr and X WD = 1,2,4,8%, from left to right, with IMF1. 
Dashed line : halo WDLF for t h = 14, 16, 18, 20 Gyr and X WD = 1.7, 8, 25, 50%, with IMF2. 
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